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ABSTRACT 



Context. It is now established that magnetic fields are present in the intra-cluster medium (ICM) of galaxy clusters, as revealed by 
observations of radio halos and radio relics and from the study of the Faraday Rotation Measures of sources located either behind or 
within clusters. Deep radio polarization observations of clusters have been performed in the last years, and the properties of the ICM 
magnetic field have been constrained in a small number of well-studied objects. 

Aims. The aim of this work is to investigate the average properties of the ICM magnetic fields, and to search for possible correlations 
with the ICM thermal properties and cluster radio emission. 

Methods. We have selected a sample of 39 massive galaxy clusters from the Highest X-ray FLUx Galaxy Cluster Sample, and used 
Northern VLA Sky Survey data to analyze the fractional polarization of radio sources out to 10 core radii from the cluster centers. 
Following Murgia et al (2004), we have investigated how different magnetic field strengths affect the observed polarized emission 
of sources lying at different projected distances from the cluster center In addition, statistical tests are performed to investigate the 
fractional polarization trends in clusters with different thermal and non-thermal properties. 

Results. We find a trend of the fractional polarization with the cluster impact parameter, with fractional polarization increasing at 
the cluster periphery and decreasing toward the cluster center Such trend can be reproduced by a magnetic field model with central 
value of few /iG. The logrank statistical test indicates that there are no differences in the depolarization trend observed in cluster with 
and without radio halo, while the same test indicates significant differences when the depolarization trend of sources in clusters with 
and without cool core are compared. The comparison between clusters with high and low temperatures does not yields significant 
differences. Although the role of the gas density should be better accounted for, these results give important indications for models 
that require a role of the ICM magnetic field to explain the presence of cool core and radio halos in galaxy clusters. 

Key words. Cluster of galaxies - Magnetic field - Depolarization — non-thermal 



1. Introduction 

Increasing attention has been devoted in the last decade 
to the presence, strength and structure of magnetic fields 
in galaxy clus ters (Ferrari etal. 2008, Carilli & Taylor 2002, 
iGovoni & Ferettil 120041 IClarkd 12004 ). In order to understand 
the physical conditions of the intra cluster medium (ICM) it 
is important to understand the properties of cluster magnetic 
fields and their interplay with the other constituents of the 
ICM. Our current knowledge on cluster magnetic fields comes 
mainly from radio observations. In some galaxy clusters mag- 
netic fields are revealed by synchrotron emission not obviously 
asso ciated with any particul ar galaxy: the so-called radio halos 
(e.g. iGiovannini et al.ll2009l) . Radio halos are wide synchrotron 
radio sources, characterized by steep spectraQ] (a >1) and low 
surface brightness (~ 1 //Jy/arcsec^ at 1.4 GHz). They are lo- 
cated at the center of galaxy clusters, and are usually found 
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to be unpola rized, the only tw o excepti ons being the c luster 
Abell 2255 (iGovoni et al.ll2005L see alsolPizzo etal.ll2010 ) and 
MACS .T0717-H3745 (iBonafede et al.ll2009ah . Sviichrotron dif- 
fuse emission is also observed at the periphery of some galaxy 
cluster (the so called rad io relics, see e.g. Bo nafede et al. 2009bl 
Ivan Weeren et al ] l2010h About 30 radio halos are known so far, 
and all of them are found in clusters with clear signature of on- 
going or recent merger a ctivity (e.g. Buote 2001, Govon let al.l 
200 ll ICassano et al.l l2010l) . Shocks and turbulence associated 
with merger events are expected to inject a considerable amount 
of energy in the ICM, that could compress and amplify the mag- 
netic field and accelerate relativistic particles, giving thus rise to 
the observ ed radio emission. We refer to Ferrari et al. (2008 ) and 
iDolag et alj (120081) for recent reviews of the subject. 
Estimates of the magnetic field strength have been obtained in 
clusters where radio halos are observed, under the minimum en- 
ergy hypothesis (which is very close to equipartition conditions) 
or by studying the In verse Compton (IC) Hard X-ray emission 
(e.g. Fusco-Femianol 12004 1^ et al.ll2009l [Aiello et al.ll2009l 
iMurgia et al.i i2010al) . Both of these methods require however 
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several assumptions on the emitting particle distribution and en- 
ergetic spectrum. In particular the assumptions required by the 
equipartition approach, in the context of radio halos, strongly af- 
fect the resulting estimates. Magnetic field values derived from 
equipartition and IC indicate that magnetic of ~ fiG level are 
spread over the cluster volum e. Recently, a more sophisticated 
approach has been applied bv Vacca et al.l (1201 0'). The authors 
investigated the magnetic field power spectrum in the cluster 
A665 by analyzing the radio halo brightness fluctuations in con- 
junction with the fractional polarization trend of radio sources at 
different impact parameters. 

Another possibility to investigate the ICM magnetic field prop- 
eities comes from the study of the Faraday rotation of sources lo- 
cated both behind and within galaxy clusters. Synchrotron radia- 
tion from radio galaxies which crosses a magneto-ionic medium 
is subject to Faraday Rotation. The direction of the polarization 
plane, ^P,,,,, is rotated by a quantity that in the case of a purely 
external Faraday screen is proportional to the square of the wave- 
length: 
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Jo 



B/fiigdl, 
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here B// is the magnetic field component along the line of sight, 
Hg is the ICM gas density and L is the distance along the line 
of sight. With the help of X-ray observations, providing infor- 
mation about the thermal gas distribution, RM studies give an 
additional set of information about the magnetic field in the 
ICM. Recent works have investigated several aspects of the mag- 
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Laing et al.l 2008) and its central strength and radial de- 



deep and multi-frequency observations of several sources lo- 
cated at different projected distances from the cluster center, and 
have thus been performed so far on a small number of clusters. 
Because of the sensitivity limits of radio telescopes, studies of a 
large number of galaxy clusters with many RM probes per clus- 
ter are still unfeasible with the current instruments. Hence in or- 
der to obtain general information on magnetic fields in galaxy 
clusters without focusing on single objects, another strategy is 
required. 

When synchrotron emission aiising from a cluster or backgro und 
source crosses the ICM, regions with similar ^F,,,,, going through 
different paths, will be subject to differential Faraday Rotation. 
If the magnetic field in the foreground screen is tangled on scales 
much smaller than the observing beam, radiation with similar 
but opposite orientation will averaged out, and the observed de- 
gree of polarization will be reduced (beam depolarization). In 
the central region of a cluster, the magnetic field is expected to 
be higher, according to both theoretical studies (see'Dolag et al.' 
2008 f or a review) and Faraday R M observation (Clarke 2004, 
Govoni et al.ll2006l iGuidetti et all 12008. Bona fede et al.ll2010l) . 
Higher value of B and ng result in higher RM, according to 
Eq. |2] The higher the Faraday RM is, the stronger the beam- 
depolarization. It is then expected that sources located in pro- 
jection close to the cluster center will show a lower fractional 
polarization c ompared to more distant ones. 
iMurgia et aP (i2004.) used the FARADAY code to simulate this 
effect, proving that a large sample of radio sources could give ad- 
ditional statistical constraints on the intra-cluster magnetic field 
properties. The advantage of this approach is that there is no 
need for multi-frequency observations. This effect can be inves- 
tigated with radio surveys at a single frequency, performed in 



full polarization mode. Such studies do not provide detailed in- 
formation on magnetic fields in specific clusters, but allow us to 
understand the average properties of magnetic fields in the ICM. 
In this paper, we present a statistical study of magnetic fields in 
galaxy clusters. Starting from a sample of clusters selected from 
die Highest X-ray FLUx Ga laxy Cluster Sample (HIFLUGCS, 
Reiprich & Bohringer 2002), we use the NRAO Northern VLA 
Sky Survey (NVSS, Condon et al. 1998) to study the polariza- 
tion properties of sources located at different impact parameters 
with respect to the cluster center 

The paper is organized as follows: In Sec.|2]we present the clus- 
ter sample, in Sec. |3] radio data are analyzed; The polarization 
propeities of the radio sources are studied in Sec. 2] and used 
to derive the average magnetic field properties (Sec. |5) with the 
help of numerical simulations. In Secs.|6]|2]and[8]the difference 
in magnetic fields in clusters with and without radio halos and 
with high and low temperature are investigated. Finally, conclu- 
sions are presented in Sec. |9] In this work we assume a ACDM 
cosmological model, characterized by Oa - 0.7, Q.,„ = 0.3 and 
Ho =70 km/s/Mpc. 



2. Selection of the clusters 

From the High est X-ray FLUx Galaxy Cluster Sample 
(HIFLUGCS) bv lReiprich & Bohringeil (l2002l) . we selected aU 
objects with L^[0.1 - 2AkeV] > 1.5 X lO'^'^erg/s. The hmit in 
X-ray luminosity is aimed at selecting massive galaxy clusters, 
and at the same time having enough clusters to build up a sta- 
tistical sample. In massive galaxy clusters the magnetic field is 
expected to be higher, so that the depolarization effect should be 
more prominent. There are 39 clusters in the HIFLUGCS sam- 
ple that satisfy this criterion. Using the Lx - M20() relation ob- 
tained by Rvkoff et al. (2008) this limit corresponds to M200 ^ 
3 X 1O'**M0//Z7o. The HIFLUGCS sample is a X-ray selected and 
X-ray flux-limited galaxy cluster sample from the ROSAT All- 
Sky-Survey catalogue. The large FOV of the ROSAT PSPC cov- 
ers most of the clusters out to rsoo, which is the radius at which 
the mean density of the cluster is 500 times that of a critical den- 
sity. The selection criterion that we used, based on the cluster 
X-ray luminosity only, guarantees that clusters in different dy- 
namical states are included in our sample. 
The HIFLUGCS sample has been used to perform several stud- 
ies regarding the scaling relations of galaxy clusters (e. g. 
IChen et al.ll2007l) . Their surface brightness profiles have been 
derived from pointed ROSAT PSPC (32 clusters) and RASS (7 
clusters) observations (Reiprich & Bohringer 2002) by using the 
standard y6-model (Cavaliere & Fusc o-Femiano 1976). From X- 
ray observations the cluster mean temperature T and the gas den- 
sity distribution have been derived. The gas density distribution 
is described by the following equation: 



ng(r) = no 1 1 -I- — 



(3) 



here r is the distance from the cluster center and is the cluster 
core radius. 

The basic properties of the c lusters in the sample, taken from 
Rei prich & Bohringeii (|2002|) and rescaled to the cosmological 
model adopted in this work, are given in Tables [1] and |2]. In the 
same Table, we list the dynamical state of the clusters as found 
in the literature, indicating clusters with merging signature (M), 
and with mild or strong cool core (CC). 
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3. Radio Data 

Radio images of the selec ted clusters were re trieved from the 
Northern VLA Sky Survey (I Condon et al.ll998l) . The NVSS was 
performed at 1.4 GHz, with a 100 MHz bandwidth, and covers 
the sky north of 6 = -40 deg. The survey was performed in full 
polarization mode, so that the total intensity radio flux (Stokes 
I) and Stokes Q and U images are provided. The images all have 
a resolution of - 45" at the FWHM and nearly uniform sen- 
sitivity. Their noise rms brightness fluctuations are cr/ ~ 0.45 
mJy/beam (Stokes I) and cry g ~ 0.30 mJy/beam in Stokes Q 
and U. Radio images for each cluster are centered on the cluster 
X-ray peak and have sizes of 20 r^x 20 re- 
The clusters SI 101, A3391, A3667, A3266, A3158 and A3112 
lying at 6 <40 deg are not in the NVSS, so that the radio sam- 
ple consists finally of 33 clusters. Coma is the nearest cluster 
of the sample, with z - 0.0232, while Abell 2163 is the most 
distant one, with z - 0.2010. Given the resolution of the NVSS 
this translates into a linear resolution ranging from 21 kpc to 
150 kpc. This guarantees that effects due to the beam depolar- 
ization should be visible at the NVSS resolution, provided that 
the magnetic field is tangled on scales smaller than these. Except 
of Abell 2163 and Abefl 2204, the clusters lie all at z<0.1, their 
mean redshift is 0.066 with a rms dispersion of 0.032. 
From Stokes U and Q images the polarization intensity and po- 
larization angle images were produced: 

P = Vt^^Tg^ 

1 U 

^ — —arctan — , (4) 

2 Q 

where P was corrected for the positive bias. 
Although the NVSS radio images have a nearly uniform sensi- 
tivity, local regions of poorer sensitivity may be present for in- 
stance because of dynamic range limitations. In order to avoid 
including in our source sample peaks of the noise fluctuations, 
only sources brighter than 5cr; in total intensity images were 
considered. We blanked the total intensity images at 5cr/, and 
used these blanked images as a mask to blank the polarization 
intensity images. We considered only sources that are wider than 
one beam area. This prevents us from including only the bright- 
est part of low brightness sources, which could bias the sample 
toward higher Fp values. We then computed the mean polarized 
intensity flux and the mean total intensity flux for each source, 
and derived the fractional polarization as Fp = P/I. 

3.1. Source detection and upper limits 

Not all of the sources detected in total intensity are also detected 
in polarization. In order to fix the flux-threshold in polarization 
to distinguish between detection and upper-limit, we proceeded 
as follows: in the NVSS images we have randomly chosen re- 
gions where no sources were detected in total intensity. These 
empty regions have been selected to have an area > than one 
beam, to mimic the same conditions of the detected sources. 
We have then computed the mean polarized intensity flux, and 
compared the distribution of "real-sources" and "noise-sources" 
(see Fig. [TJ. The threshold was fixed at the value of P where 
the number of "noise-sources" is less than 10% than the num- 
ber of "real-sources". This ensures that we are safely dominated 
by detections at every R The threshold corresponds to P, -0.45 
mJy/beam so when P < P, mJy/beam a 1-cr upper-limit for Fp 

was derived: Fp < Jjj^^. 
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Fig. 1. Number of sources and noise as a function of the observed 
polarized intensity. 

The final sample of sources consists of 696 sources, 334 of 
which are upper Umits. 

4. Depolarization Analysis 

The large sample of sources obtained by the NVSS catalogue al- 
lows us to investigate how the depolarization changes as a func- 
tion of the projected distance from the cluster center 
A radial decline of the magnetic field proportional to n'' is ex- 
pected from magneto- hydrodynamical simulations performed 
with diff'erent codes (Dolag et al.' '2005' 'Dolag et al.' '19991 
Bru2genetal. l2005l [Dubois & Teyssier 2008 , ' Collins et aT| 
2010) as a result of the compression of thermal plasma during 
the cluster gravitational collapse. In addition , the comparison be- 
tween thermal and radio profile brightness (iGovoni et al 1 1200 Ih 
also suggests that the non-thermal component follows the distri- 
bution of the thermal gas. Here we will assume that the magnetic 
field radial profile is proportional to the gas density profile, ac- 
cording to: 

B(r) = <Bo)(^j', (5) 

where 77 is a parameter Given this assumption, the different val- 
ues of rc for the clusters in our sample have to be accounted for 
The gas density profile (Eq. |3]l is almost constant for r < and 
decrease faster as r > r^. We have normalized the distance of the 
sources to the cluster core radius, defining r„o„„ - r/rc. In the 
upper panel of Fig.|2]the depolarization percentage as a function 
of r„orm is reported for all the sources and the upper limits in our 
sample. 

4.1. Fractional polarization trend 

The sample of sources that we have obtained is affected by 
the presence of non-detections, or censored data. Censored data 
points are those whose measured properties are not known pre- 
cisely, but are known to lie above or below some limiting sensi- 
tivity. Here a non-detection still gives us the useful information 
that the P flux (and so the fractional polarization) is below a cer- 
tain threshold. We have thus a left-censored sample of data. Our 
aim is to derive the depolarization trend of the sample of sources 
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we have selected, taking care of the information given by cen- 
sored data. 

An extensive field of statistics called "survival analysis" of "life- 
time data" exists to address problems of this kind. A very useful 
statistical estimator in the survival analysis is the Kaplan Meie r 
(KM) estimator, that is described bv lFeigelson & NelsonI (Il985h . 
It is a non parametric maximum likelihood type estimator of the 
distribution function F(t), with t being a generic variable, in our 
case Fp. It is usually expressed in term of the survival function 
S(t). Let xi < X2 < .... < X, denote the distinct, ordered, ob- 
served values, the KM estimator is given b>0: 

S KM(t) = ProbiT > t) = 1 - Fit) (6) 
= ri/..^. <, ^1 ~ ~j ' ■^hen t > x\ (7) 
- 1, when t < x\ (8) 

with rii being the number of objects (detected or undetected) 
> X;, di are the number of objects at value x,, and 5, = 1,0 if 
Xi is detected or undetected respectively. 

Using the KM estimator, we have derived the distribution func- 
tion of our data set in each bin. In the lower panel of Fig. |2]the 
median of the distribution function in each bin is reported. The 
error bars refer instead to the 16th and 84th percentile of the dis- 
tribution (those that include 68% of the data in the distribution). 
A trend is detected in the Fp going from the cluster center to the 
cluster outskirts (Fig.|2]i. The value of Fp in the first bin, i.e. for 
rnorm <1 is 0.0046, and increases gradually in the outer bins. The 
effect of the cluster on the observed Fp is clear out to ~5 core 
radii, while in the outer bins, the mean Fp is scattered around 
a constant value. It reaches values going from 0.047 to 0.058 in 
the 7th to 10th bin. As the distance from the cluster center in- 
creases, RM decreases and the observed Fp is no longer affected 
by the presence of the ICM. 

The value that we find in the outer bin is interpreted as due to de- 
polarization intrinsic to the structure of the radio source as seen 
by the radio-telescope at this resolution and frequency. We as- 
sume that there is no bias in the intrinsic depolarization in all 
the sources in our sample. The intrinsic fractional polarization is 
assumed to be the same for sources regardless of their projected 
distance from the cluster center, as found e.g. in the Coma clus- 
ter (Bonafede etal. 2010). 

It follows that internal depolarization should not affect the ob- 
served trend of Fp versus r„o„n, but rather act like a constant 
normalization factor 

The observed trend indicates that magnetic fields are common 
constitu ents of gal axy clusters, in agreement w ith the results by 
IClarkd (§004) and iJohnston-HolUtt et al] (12004) . who analyzed 
the Faraday Rotation Measures of sources located behind and 
within clusters. 

5. Magnetic fields properties 

The observed trend of the fractional polarization Fp versus the 
projected distance from the cluster center can be used to con- 
strain the magnetic field properties. Both the strength and the 
morphology (/. e. the power spectrum) of the magnetic fields af- 
fect the Faraday RM of radio sources, and the consequent beam 
depolarization that we are investigating here. 

^ Here the KM estimator is expressed in the case of right-censored 
sample. To obtain an estimator of F{i) for left-censored samples we 
have computed P(T < t) = S km(M - t)), with M being the maximum 
of ti. 



We use 3D simulation s of random magne tic fields, performed 
with FARADAY code ("Murgia e t ani2004 to simulate the ex- 
pected trend of Fp for different magnetic field configurations. In 
this Sec. the magnetic field modeling is briefly presented, and 
the comparison with the observed sample is performed. 

5. 1 . Magnetic field modeling 

The FARADAY code dMurgia et al.ll2004l) simulates 3D models 
for the magnetic field in galaxy clusters, and computes the ex- 
pected RM for radio sources within/behind the cluster. 
Numerical simulations start by considering a power-law power 
spectrum for the vector potential A in the Fourier domain: 

\Ak\^ (X (9) 

and extract random values of its amplitude A and phase (p. The 
amplitude A is randomly extracted from a Rayleigh distribution 
(in order to obtain a Gaussian distribution for the real magnetic 
field components), while <p varies randomly from to 27r. Here 
k is the wave vector in the Fourier domain. The correspondent 
quantity in the real space is A = ^. The magnetic field compo- 
nents in the Fourier space are obtained by computing the cross- 
product: 

B(k) = ikxA(k). (10) 

Finally, the field components B, in the real space are derived us- 
ing a 3D Fast Fourier Transform inversion. The numerical mag- 
netic field model is then divergence-free, Gaussian random, and 
isotropic. The power spectrum normalization is set such that the 
average magnetic field strength follows Eq. |5] This operation 
is performed in the real space domain, to reduce the computa- 
tional burden. The same approach has already been used in sev- 
eral works (Murgia et al. 2004, Laing et al. 2008, Guidett TetaD 
20081 iBonafede et al.ll201QllVacca et alJl2010() . 
The adopted power spectrum power law introduces three free pa- 
rameters: «, A„„„ and A,„a_x- Two more free parameters are intro- 
duced once the magnetic field profile (Eq.|5]l is considered: (Bq) 
and T]. Two main degeneracies affect these parameters: the first 
one concerns the power spectrum, and in particular the values of 
n and A„„„, and another degeneracy affects (Bq) and rj. Different 
combinations of n and A„ax as well as (Bo) and 77 lead to similar 
values of fractional polarization, since it depends mainly on the 
amount of Faraday Rotation Measures that originates in the ICM 
causing the beam depolarization. 

Given the number of free parameters, we have fixed the magnetic 
field radial profile slope (7/), and for three different magnetic field 
power spectra we have investigated the magnetic field strength 
at the cluster center 

The magnetic field power spectrum is assumed to be 
Kolmogorov-like. Different independent analyses of RM data 
(Vogt & EnBlin 2005, Guidetti et al. 2008, Bonafede et al. 201§ 
indicate that the magnetic field power spectrum observed in 
galaxy clusters is in good agreement with the power law ex- 
pected by the Kolmogorov theory, which in our 3D notation 
corresponds to n = 11/3. We note that the theory developed 
by Kolmogorov is applicable incompressible un-magnetized and 
uniform fluids, so that its application to the case of the ICM 
of galaxy clusters is all but obvious. Nonetheless, observational 
data and cosmological simulations indicate that it is a good de- 
scription of the pseudo-pressure fluctuations (Schueck er et al.l 
2004), velocity field ( Vazza^t al. 2009) a nd magnetic field ( 
:yogL& EnBlin 2005, G uidetti et a l. 2008, Bonafede et al.ll2010h 
in the ICM. Although other spectral slope are not excluded by 
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Fig. 2. Upper panel: Fractional polarization for the sources belonging to the cluster sample as a function of the projected distance 
from the cluster center normalized by the core radius. The arrows indicate upper limits. Lower panel: the median of the KM estimator 
in every bin is shown, error bars refer to the 16th and 84th percentile of the KM distribution. 



data, we will adopt here a Kolmogorov like power spectrum. The 
power spectrum maximum and minimal scale were set to repro- 
duce a magnetic field power spectrum having an auto-correlation 
length Ab ~164 kpc (model 1), ~ 25 kpc (model 2) and 
Afi ~ 8 kpc (model 3). We then adopted A„„„ =8 kpc, much 
smaller than the beam linear size. The values of A,„ax are then 
512 kpc (model 1), 64 kpc (model 2), and 14 kpc (model 3). The 
value of T] is set to 0.5. This choice is motivated by recent works 
(iGuidetti et al.ll2008l iBonafede et"an 2010) that found the best 
agreement between simulations and observations when the mag- 
netic field energy density follows the thermal gas density ( i. e. 
77=0.5). We then investigated different values of {Bo}=l, 5, and 
lOfiG. 



5.2. Comparison with observed data 

Numerical simulations of the magnetic fields have been per- 
formed in a cubical box of 1024^* pixels, with a pixels-size of 2 
kpc. The simulated magnetic field is periodic at the grid bound- 
aries, and the computational grid has been replicated to obtain 
a simulated field of view of ~ 20 (rc)x 20 (re), with (re) being 
the mean value of for the clusters in the sample weighted by 
the numbe r of sources in each cluster (~ 270 kpc). 
Following iMurgia et al.l (|2004|) our aim is to compute, for dif- 
ferent magnetic field models, the ratio between the polarization 
obtained after the rotation of the plane of polarization and the 
intrinsic polarization at a given frequency. This is the quantity 
that can be compared with the observed Fp, once the effects in- 
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troduced by observations are considered. 

We simulated different magnetic field models and derived the 
expected Fp for sources located halfway through the cluster at 
increasing distance from the cluster center We proceeded as fol- 
lows: 

- We simulated different magnetic field models, with (Bo) = 
1, 5 and 10 yuG, and = 8, 25, and 164 kpc. 

- The intrinsic polarization f of the sources is assumed to be 
spatially uniform and ordered (^F^, constant) and its value is 
set to reproduce the observed value of Fp at r„o™ =10. 

- According to Eq. |2] we derived the synthetic RM images 
expected for the different magnetic strengths. The integra- 
tion path goes from (cluster center) to 10 {re). We used 
no - 0.004 cm"^, fi - 0.67 that are the mean values of our 
sample weighted by the number of sources that are in each 
cluster We produced simulated RM images on a 5400x5400 
kpc field of view, with a resolution of ~2 kpc. 

- We then computed the expected polarization plane direction 
according to Eq. 1 and taking into account the effect due to 
the NVSS finite bandwidth. 

- Pint and "^obs were then converted into the Stokes parameters 
Q{A) and U(X) following Eq.g] 

- To the simulated images of Qi/i) and U{A) we added a 
Gaussian noise having rms=0.3 mJy/beam, as the mean 
NVSS noise for Q and U images. Q and U images were then 
convolved with a Gaussian beam having major and minor 
axis ^ 42 kpc. This is the linear size corresponding to 45" 
at z=0.048 that is the median redshift of the cluster sample 
weighted by the number of sources in each cluster 

- Finally, the synthetic Q and U images were transformed back 
to P using Eq. |4] 

- The trend of simulated and observed P versus rnom have 
been compared. 

The comparison between the simulated and observed Fp trend 
obtained by averaging data from the clusters in the sample is 
performed under the following assumptions: 

1 . the magnetic field is assumed to be a Gaussian random field. 
This can be incorrect if the magnetic field is intermittent, 
and characterized by filaments. However, this is a common 
assumption in this field, and is mainly motivated by the fact 
that current data are gener ally too sparse to d etermine higher 
order correlation (see e.g. lLaing et al.ll2008l) . 

2. The magnetic field phases are assumed to be random, 
which means that the magnetic field is assumed to have no 
preferred orientation. Again, this may be no t the case when 
individual clusters are analyzed ( e. g. M 84. lLain^ & Bridlj 
[1987) but we are interested in the magnetic field properties 
averaged over large volumes of several clusters. 

3. The power spectrum fluctuations are assumed to vary with 
the thermal gas density in the ICM, according to Eq. |5] 



Under these assumptions, the model that best reproduces the ob- 
served trend of Fp, among those considered here, has (Bq) - 

The simulated Fp trends for the different magnetic field models 
are shown in Fig.|3] We are aware that we are probably averaging 
the Fp produced by slightly different magnetic field models and 
configurations, and that more refined and cluster specific models 
should be considered if precise values are needed. Nonetheless, 



the best-fit model gives a reasonable order of magnitude estimate 
for the average magnetic field properties in galaxy clusters. 

In our best fit model with a central magnetic field (Bo) = 
5//G, the mean magnetic field over the central Mpc cube, B^iMpc), 
resulting from our best fit model is ~2.6 juG. A more detailed 
analysis has been performed for some of the clust ers in our sam- 
ple: Abell 2255 (Govo niet al1l2006h . Hydra- A (Vogt & Enfilin| 
2005, Laingetal. 200l) and Coma ffionafede et al. 201(17 
Given the degeneracy between (Bo) and we compare the mean 
magnetic field over the central Mpc cube, found by the above 
mentioned autho rs with the one obta ined here. In A2255 and in 
the Coma cluster lGovoni et al.l (|2006|) and lBonafede et al.l (l20I0h 
find B( 

I Mpc) •^ 1.2 //G, a nd B(\mpc) ^2 jjG respectively, while 
the analysis bv lVogt & E nfilin (2005) yields B(\Mpc) ~1 A'G for 
Hydra A. They are compatible with the value obtained here. 

6. Dependence of magnetic field properties in 
clusters with and without radio halo 

The sample of clusters that we have selected comprises both viri- 
alized and merging systems (see Tab. |2]and references therein). 
In a fraction of clusters, with signatures of recen t or on-going 
merg ers, a radio halo h as been detected (see e.g. IVenturi et al.l 
20081 iGiovannini et al ][2009). The origin of these radio sources 
is still un-known, and two main classes of models have been pro- 
posed in the literature: 

- primary or re-acceleration models: in which electrons are re- 
accelerated in situ through second-order Fermi mechanism 
by ICM turbulence developing duri ng cluster mergers (e.g. 
iBrunetti et al.ll200lHPetrosianll200ll) : 

- secondary or hadronic models: in which electrons originate 
from hadronic collisions between the long-living relativistic 
protons in the ICM and thermal ions (e.g. JDennison 1980). 

One possibility to explain the presence of radio halos in only a 
fraction of clusters is to assume different magnetic field strength 
in the ICM of clusters with and without radio emission. The 
magnetic field could be amplified during cluster mergers, as 
indicated by MHD cosmological simulations (see review by 
iDolag et alJf2008) . and then, it could potentially be dissipated 
once the cluster turns back to the equilibrium state. 



6.1. Clusters with and without radio halo 

With the aim of investigating a possible difference between the 
magnetic field properties in clusters with and without radio ha- 
los, we have divided our initial sample in two sub-samples, that 
contain clusters with radio halos (labeled with H in Tab. |2]i and 
clusters where a giant radio halo (i.e. with size ~1 Mpc) has not 
been detected so far It is worth noting that no upper limits have 
been put on the radio emission of clusters where a radio halo has 
not been detected, so that additional care is required to avoid, or 
at least limit, obvious observational bias. 

Most of radio halos have been discovered with follow-up ob- 
servations of NVSS candidates. Since the NVSS cannot detect 
large (Mpc) scale emission in clusters at z<0.044, we exclude 
from our initial sample the clusters with such a redshift. Note 
that this is a conservative approach, since other single-dish in- 
struments would be clearly able to detect radio emission in these 
very nearby objects, as in the case of the Coma cluster 
In addition, the radio power of the observed radio halos is 
usually found to correlate with the cluster X-ray luminosity 
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Fig. 5. Fractional polarization versus the projected distance from the cluster center, normalized by the core radius, for sources in 
clusters that host a radio halo (left), and that do not host a radio halo (right). Arrows indicate upper limits. Bottom panel: crosses 
refer to the median of the KM estimator is each bin, bars indicate the 14th and 84th percentile. Due to the low number of points in 
each bin the width of each bin is adaptively set so that at leat 3 detections fall within a bin. 



(e.g. lLiang et alJl2000llFeretti et al.ll2000l) . Our sample includes 
clusters with X-ray luminosities going from 1.5 to 20.8 xlO'*'* 
erg/s. Clusters with the lower X-ray luminosity could host ha- 
los of lower radio power, that are undetected by present sur- 
veys (see e. g. IClarkg|2005h . It is worth mentioning that recently 
two radio halos have been fo und in clusters wit h low X-ray 
luminosity, namely 021 7-H070 (iBrown et al.ll201 ih . and A1213 
dGiovannini et alj |2009), whose radio power is an order of mag- 
nitude or more above the extrapolation of the Lx-Pr correlation 
at low X-ray luminosities. However, we decided to adopt a con- 
servative approach, and to compare clusters in the same X-ray 
luminosity range. In Fig.|4] the X-ray luminosity of the clusters 
in the energy range 0.1- 2.4 keV is shown. Clusters are repre- 
sented with different colors and symbols, depending on the pres- 
ence of radio halo and on their dynamical state. Objects with 
z<0.044 have been removed. Clusters with and without radio 
halo populate the Lx range uniformly. 

Among the clusters that are known to host a radio halo, 
the cluster A3562 is the one with the lowest X-ra y luminosity: 
1.67-tl O"^"^ erg/s. The radio halo was discovered bv ' Venturi et al.l 
JlOOOh by inspecting NVSS images. We will then consider in the 
following comparison only clusters with Lx higher than the one 
of A3562. 

Although these two additional criteria should prevent us from 
obvious observational bias, it is worth mentioning that radio ha- 
los with very steep radio spectra and/or low radio brightness 
could still be present in the cluster labeled as "radio-quiet". Only 
future instruments, such as LOFAR, will be able to provide a 



definitive answer Nonetheless, it is interesting to investigate the 
properties of the ICM magnetic fields in these two sub-samples, 
according to our present knowledge about their radio emission. 

The fractional polarization for the two sub-samples (radio- 
halo and radio-quiet) is shown in Fig. |5] The halo sample con- 
sists of 374 sources (194 of which are upper limits), while the 
radio-quiet sample consists of 243 sources (112 upper limits). 
We want to test the null hypothesis that the two populations, sub- 
ject to left-censoring, have the same distribution. We performed 
the logrank statistical test, that is a non-parametric test, widely 
used in astronomy when censored populations have to be com- 
pared. 

According to the logrank test, let Tij being the observed value 
from distribution / - 1,2, with /' = l...A^,-. Let then yi < y2 < 

< yr with r < N = Ni + N2 denote the uncensored 

values in the combined samples. Given two independent ran- 
dom samples drawn from populations having distribution func- 
tion^ Fi(t) - P(Tij > t) we want to test the null hypotheses 
Ho : Fi(t) = ^2(0^ for all t. Let dij denote the number of objects 
from sample / = yj, dj the number of objects from both samples 
- yj, tii j the number of objects from sample / > yj and similarly 
rij the number of objects from both samples > yj. The quantity 



djn\ 
nj 



(11) 



^ Here the logrank test is expressed in the case of right-censored sam- 
ple. 
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Fig. 6. Top panel: Fractional polarization versus the projected distance from the cluster center, normalized by the core radius, for 
sources in cool core clusters (left) and non cool core clusters (right). Arrows indicate upper limits. Bottom panel: crosses refer to 
the median of the KM estimator is each bin, bars indicate the 14th and 84th percentile. Due to the low number of points in each bin 
the width of each bin is adaptively set so that at leat 3 detections fall within a bin. 



for large A^, under Hq, is approximately normally distributed 
with zero mean and variance cr^. Hence, Hq is rejected at level a 
if 



where Za/2 is the score such that the area under a standard nor- 
mal curve over the interval [-Za/2,Za/2] is l-a. We refer again to 
iFeigelson &Nelsonl(ll985h for an illustrative explanation of the 
test with application to astronomical data. Our null hypothesis 
is that the two populations of Fp in clusters with and without 
radio halo are different realization of the same sample. The null 
hypothesis is accepted with high significance: P = 0.99. We can 
conclude that the depolarization trend for sources seen through 
clusters is the same in these two sub-samples. Since the amount 
of depolarization depends on the Faraday rotation in the external 
screen, there are two different possibilities to explain our results: 

1. Both magnetic field and gas density distributions in cluster 
with and without radio halos have similar properties, or 

2. magnetic field and gas density are both different in clusters 
with and without radio halo, but conspire to give rise to the 
same Fp. (i.e. higher magnetic field and lower gas density in 
one sample and lower magnetic field and higher gas density 
in the other one) 

The mean values of no, weighted by the number of sources, for 
the radio halo and radio quiet samples are 0.003+0.001 cm"-* 
and 0.005±0.010 cm"^ respectively. Since they are compatible 



within errors, we can conclude that there is no evidence for a 
different magnetic field in clusters with and without radio halos. 

6.2. Analysis of the possible contamination by merging 
effects 

Processes related to merger events could change the magnetic 
field structure in galaxy clusters, increasing the field auto- 
correlation length and thus also the Fp of the observed sources. 
Also, if merger events amplify the magnetic field strength, a 
lower value of Fp is expected. These two effects could conspire 
to give the same Fp distribution in the radio halo and radio quiet 
samples even if magnetic field is higher in the radio halo sam- 
ple, because of the larger auto-correlation length of the magnetic 
field. In order to prevent this, we have repeated the analysis de- 
scribed in Sec.|6]considering only merging systems in the radio 
quiet sample. The radio quiet merging sample consists of 203 
sources (94 upper limits), and the logrank test yields a probabil- 
ity of 0.70 that the radio halo and radio quiet merging sample are 
drawn to the same intrinsic population. We can conclude that the 
depolarization trend is the same in merging clusters, regardless 
of the presence of a radio halo or not. 

6.3. Magnetic field and radio halo origin 

The results we have obtained in the previous sections can also 

be used as a test for the radio halo formation models. 

iBrunetti et al.l (l2009i) have analyzed the well-known correlations 
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Fig. 7. Fractional polarization versus the projected distance from the cluster center, normalized by the core radius, for sources in 
clusters with T> 7 keV (left), and clusters with T< 7 keV (right). Arrows indicate upper limits. Bottom panel: crosses refer to the 
median of the KM estimator is each bin, bars indicate the 14th and 84th percentile. Due to the low number of points in each bin the 
width of each bin is adaptively set so that at leat 3 detections fall within a bin. 



between Radio halo power at 1.4 GHz, P| 4 and the cluster 
X-ra y luminosity, Lx (e. g. iLiang et al.ll200"oi iGiovannini et al.l 
120091) in a sample of X-ray selected clusters. Among this sam- 
ple only =s I of clusters host radio halos, while for the other 
clusters o f the sample upper l imits on their radio emission have 
been put dBrunetti et al.ll20() 7'). It must be noted that the sample 
of clusters that they have analyzed is composed of both merging 
and dynamically relaxed clusters, and that the upper-limits have 
been computed assu ming a value for the radio spectral index. 
iBrunetti et al.l (l2009l) conclude that in the hadronic scenario the 
energy density of the magnetic field in clusters without a radio 
halo should be at least 10 times lower than in those that host ra- 
dio halo emission. 

More recently, iKeshet & Loebl (1201 Oh have proposed a model 
within the hadronic scenario, in which cosmic rays diffuse away 
from their sources (supposed to be supernovae), whereas the 
magnetic fields are amplified by mergers in clusters with radio 
halos. The observed bi-modality in the Pi^-Lx plane is then at- 
tributed to a bi-modality in magnetic field strength for clusters 
with and without radio halos. Radio halos are associated with 
clusters having strong B, i. e. B » Bcmb in the radio emitting 
region, whereas clusters with B « Bcmb all over their volume 
would be radio quiefl 

The result obtained in the previous sections indicates that the 
bi-modality in the Pi^-Lx plane cannot be attributed to a bi- 
modality in magnetic field properties. 



The result we find here is based on a statistical approach, how- 
ever, previous works have already investigated the magnetic field 
properties in clusters either with and without radio halo, find- 
ing results consistent with those we obtain here. For example, 
magnetic fields of seve ral nG were inferre d in clusters with- 
out ra dio halo by e.g. 'Clarke et al.' ("2001"), and iMurgia et"an 
(l2004l) . and recently, Govoni et al. (2010) found that clusters 
seem to follow a common Sx - RM distribution, Sx being the 
X-ray surface brightness, independently on their radio proper- 
ties. However, this is the first time that a statistical comparison 
between magnetic fields in clusters with and without radio emis- 
sion has been performed and our result is a further indication 
that hadronic models are difficult to recoricile wi t h present data 
(IBrunetti et al] l2008l and iDallacasa et al.1 120091; iDonnert et all 
iJeltema & Profumol2010l) ! 



* Bcmb = 3.2(1 + zf^iG is the equivalent magnetic field strength of 



the cosmic microwave background (CMB). 



7. Dependence of magnetic field parameters on tlie 
presence of a cool core 

Galaxy clusters are often divided, according to their X-ray emis- 
sion, into cool core and non cool core clusters (e.g. Hudson et aT| 
20101) . The former class is characterized by a bright central 
peak in the cluster X-ray surface brightness, a high central gas 
density, and a positive gradient of the metal abundance profile. 
The core is characterized by a radially increasing temperature. 
Cool core clusters are believed to be systems in dynamical 
equilibrium, while the energy released by cluster mergers could 
transform cool core clusters into non cool core systems (see 
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Fig. 3. Simulated (lines) and observed (crosses) depolarization 
trends for models with different =164,25 and 8 kpc (top, 
middle and bottom panel respectively). Different colors refer to 
different magnetic field strengths, as reported in the figures. 



e.g. lMotletaljr2004 iRossetti & MolendillloTol and references 
therein). In the core, the cooling time of the gas is much shorter 
than a Hubble time, but the lack of strong cooling lines revealed 
by X-ray spectroscopy of Chandra and XMM-Newton, indicates 
that the cooling of the gas must be inhibited. Several processes 
could in principle be responsible for that, such as the presence 
of heating sources (e.g. iBinne v & Tabor 1995), or therma l 
conduction from larger radii, e.g. Naravan & Medvedevll200ll) . 



Fig. 4. X-ray luminosity in the 0. 1-2.4 keV band for the clusters 
in the sample here analyzed. Different symbols and colors re- 
fer to different sub-samples. Clusters that host a radio halo (red 
diamonds); clusters where a radio halo has not been detected so 
far, but with signature of on-going merger according to the litera- 
ture (cyan triangles), other clusters (squares). The cluster A1650, 
whose merging state is debated, is signed with a blue asterisk. 



Recently, several authors have investigated the role that mag- 
netic fields could have in this context (e.g. iRuszkowski & OhI 
[2010, Parrish et al. 2010). 

We will compare here the depola rization trend of sources in 
clusters with and without cool core. IChen et al] (l2007i) analyzed 
the basic properties of the clusters in the HIFLUGCS sample 
(see Tables 1 and 2 of that paper) and divided the objects in 
moderate cool core, pronounced cool core and non cool core. 
Both moderate and pronounced CC are labeled as CC in our 
Tab. HI 

We have divided our initial sample into two sub-samples. 
Since cool core clusters are characterized by the presence of 
a cD galaxy sitting in the cluster center, while the same is not 
always true for non cool core cluster clusters, we considered 
in the non cool core sample only clusters that have a radio 
galaxy at their center (A1736, A2163, A3558, A3562, COMA, 
and ZWCL1215). The cool core sample consists of 66 radio 
sources (30 of them are upper limits), while the non cool 
core samples comprises 295 sources (165 upper limits). Fp 
data are shown in Fig. |6] We want to test the null hypothesis 
that the two populations, subject to left-censoring, have the 
same distribution. We performed the logrank statistical test. 
The result of the logrank tests applied to the samples of cool 
core and non cool core clusters is that the null hypothesis has 
a small significance (P=0.15). Although the null hypothesis 
cannot be rejected at 95% confidence level (corresponding to 
P=0.05), that is the threshold commonly assumed in statistics, 
the low value of P indicates that the two samples are likely to be 
realizations of intrinsically different populations. This could be 
due to a different magnetic field or to the different gas density 
that characterize the sample (we find a weighted mean of 
no =0.017±0.013 for the cool core sample and no =0.003±0.001 
for the non cool core sample), or a combination of the two. 
Increasing attention has been devoted in recent years to the 
role that different magnetic field configurations could have 
in explaining the observed bi-modality among cool core and 
non-cool core clusters (see e.g.J Parrish et al. 2010, Parrish et 
120091 IRuszkowski et all l2010l and .Ruszkowski & Oh ,201(1 
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The main idea is that magneto-thermal instabilities could lead 
to different magnetic field configuration in the ICM depending 
on the thermal properties of clusters. When VT ■ g < 0, i.e. 
when the gravity field g is anti-parallel to the temperature 
gradient in the ICM, heat-driven flux buoyancy instabilities 
(HBI) are developed. The HBI reorient the magnetic field 
lines to be perpendicular to the temperature gradient, thus 
reducing the effective conductivity of the plas ma and preventing 
the coohng catastrophe dRuszkows ki & Oh | |2010). Recently, 
iRuszkowski & OhI (120101) and P arrish et all (12010) have shown 
that modest levels of turbulence (~ 100 km/s) can suppress 
the HBI, resulting in a quasi-stable thermal equilibrium, with 
isotropically tangled magnetic field lines. However lower levels 
of turbulence mixing are insufficient to suppress HBI, resulting 
in a thermal runaway and leading to a cool core cluster A 
different magnetic field configuration would then be expected in 
clusters with and without cool cores, and the combined effects 
of HBI and turbulence would explain how minor and major 
mergers could disrupt the cool core. 

At the aim of quantify how much of the observed Fp is due to 
the gas density and magnetic field inside the core, we integrated 
numerically Eq. |2] from to one core radius and from 1 to 
10 core radii, assuming a magnetic field with {Bq) =5/vG, 
«o =0.017 cm"^, /3 - 0.6, =90 kpc, that are the weighted 
means for sources in the cool core sub-sample. The ratio of the 
two quantities gives 0.79, indicating that the RM resulting from 
the core is a consistent fraction of the total amount which is 
responsible for the observed Fp. 

The results we have obtained in this section is not conclusive, 
but suggest a possible difference in the magnetic field properties 
of clusters with an d without cool cores , thus s uppo rting the sce- 
nario proposed by IRuszkowski & OhI (1201 Ol) and iParrish et al.l 
( 1201 0). Unfortunately, these data are not suitable for further in- 
vestigation (i.e. to disentangle the effects of the gas density and 
to better investigate if and how the difference in the magnetic 
field properties of these two samples is due to different magnetic 
field strengths and structure) but indicate that a more detailed 
analysis of the magnetic field in these two samples could offer 
important information about the interplay of magnetic fields, 
thermal conduction and the bi-modality cool core, non cool core 
clusters. 



8. Magnetic field strength and cluster temperature 

Recent works have analyzed a possible connection between 
the magnetic field strength and the cluster gas density and 
temperatur e. These studies a r e based on cosm ological simula- 
tions (e.g. |DolagiitaD |2Q^ 

l2005h or plasma physical considerations (iKunz et al.l2010h . 
SPH simulations predict that the mean magnetic field strength 
varies with the ICM mean temperatur e according to^B oc T^. 
A shallower trend has been found by, iKunz et aP (1201 Oh who 
found B oc T^/"^. 

Here we investigate a possible connection between the observed 
Fp and the cluster temperature. The selection criterion we 
have used, based on high X-ray luminosity, naturally favors the 
selection of hot clusters. Nonetheless, clusters in our sample 
span a good range of temperatures, going from 3 to 13 keV. We 
have divided our initial sample in two sub-samples having T > 
7 keV and T <7 keV respectively. The cut in temperature is set 
to have almost the same number of sources in both samples. 
The high and low T sub-samples consist of 305 sources (163 
upper limits) and 391 sources (171 upper limits), respectively. 



The trend of Fp versus r„o„„ for these two sub-samples is shown 
in Fig.|2l The logrank test yields that the null hypotheses should 
be accepted with P=0.64, thus indicating a common origin for 
the two population. The mean central gas density no for the two 
samples is 0.003+0.010 (T>7 keV sample) and 0.004±0.015 
{T < 1 keV sample), meaning that the logrank test on Fp is 
actually a test on the magnetic field properties, while different 
values of the gas densi t y sho uld not play an important role. 
Recentlv lGovoni et al.l ( 12010 ) have analyzed the Faraday RM of 
a sample of sources that belong to hot nearby galaxy clusters, 
and used these data together with literature ones in order to 
investigate a possible connection between the magnetic field 
strength and the cluster mean temperature. Clusters are divided 
into three samples having T<4 keV, 4<T<8 keV and T>8 keV. 
The data analyzed do not show evidence for such a correlation, 
indicating that a possible connection between the magnetic 
field strength and th e gas temperature, if present, is very weak 
dGovoni et al.l l2010l) . As previously noted, the temperature 
threshold of 7 keV has been chosen arbitrarily to divide clusters 
in to two s ub-samples. In order to compare our results with those 
of lGovon i et al. (2010), we repeated the test comparing clusters 
with T<4 and T >8 keV. There are 75 sources (28 upper limits) 
in clusters with T<4 keV and 305 sources (163 upper limits) 
in clusters with T>8 keV. The logrank test indicates that the 
null hypothesis can be accepted with P=0.66. These data do not 
indicate a possible difference in the magnetic field properties of 
clusters depending on their temperature, but due to the small 
number of sources in the T<4 keV sample, a more detailed 
analysis would be required. We note however that these results 
are in agreement with those obtained by RM data (Govoni et al. 
2010). 



9. Conclusions 

We performed a statistical analysis of the fractional polarization 
of radio sources in a sample of X-ray luminous galaxy clus- 
ters with the aim of studying the properties of the intra-cluster 
medium magnetic field. We used data from the NVSS to search 
for a trend in the fractional polarization with distan ce from the 
cluster center, following the approach proposed by Murgia et alj 
(l2004l) . The sample of clusters we have selected comprises both 
cool core clusters and merging systems. We have detected a clear 
trend of the fractional polarization, being smaller for sources 
close to the cluster center and increasing with increasing dis- 
tance from the cluster centers. The low fractional polarization 
in sources closer to the center is interpreted as result of higher 
beam-depolarization, occurring because of the higher magnetic 
field and gas density in these regions. This result confirms that 
magnetic fields are ub iquit ous in galaxy clusters , as al ready 
found by IClarkd (|2004 and iJohnston-Holhtt et al.1 (|2004 . Our 
results can be summarized as follows: 

- We used 3D magnetic field simulations performed with the 
FARADAY code to search for the magnetic field model that 
best reproduces the trend of the fractional polarization with 
the distance from the cluster center We assumed that the 
magnetic field power spectrum is Kolmogorov-like and that 
the magnetic field energy density decreases with radius as 
the gas thermal energy density. Under these assumptions, we 
found that a magnetic field with a central value of 5/iG gives 
the best agreement with the observed fractional polarization 
trend. 
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- We investigated possible differences between the Fp trend 
observed in clusters with and without radio halos. The 
logrank statistical test indicates that the two sample of 
sources very likely belong to the same intrinsic population 
(P=0.99). Magnetic fields in galaxy clusters are then likely 
to share the same properties regardless of the presence of 
radio emission from the ICM. This result poses problems 
for the "hadronic-models" for the origin of radio halos, that 
requires a difference in magnetic field strengths in clusters 
with and without radio halos. 

- We searched for possible differences in the magnetic field 
properties in clusters with and without cool core. The 
logrank test indicates that the Fp distribution observed 
in clusters with and without cool cores is likely to be 
different (the null hypothesis of the two samples belonging 
to the same population has a low significance: P=0.15). 
This is expected by recently proposed models that explain 
the cool core and non cool core bimodality in as due to 
different magnetic field configurations in clusters with and 
without cool core (Parrish et al. 2009, Parrish et al. 2010, 
lRuszkowskietalJ l2010. and Ruszkow ski & OhI 120100 . The 
results obtained here must be taken with cautions since the 
role of the different gas densities in the two samples is not 
easy to quantify, and could play a crucial role. Deeper radio 
observations would be required to properly test these models. 

- We searched for a possible dependence of magnetic fields 
from the cluster mean temperature. Clusters were initially 
divided into two samples having T > 7 keV and < 7 keV. 
The result of the logrank test yields that the null hypothesis 
that the two samples are intrinsically drawn from the same 
population is accepted with P=0.64. In order to compare the 
Fp analysis with the recent result by Govoni et al. (2010), 
we repeated the test by dividing the clusters into two sam- 
ples having T> 8 keV and T<4 keV. The logrank test indi- 
cates now that null hypothesis can be accepted with P=0.66. 
Due to the smaller number of sources in the T<4 sample 
this result would need further investigation, but indicate that 
clusters with different temperatures share the same magnetic 
field properties. 
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Appendix A: Investigating possible caveats 

We will investigate here possible caveats that could affect our 
analysis of the fractional polarization trend (Sec.@}. 



"wide-angle tail" type sources are present in the cluster. The lo- 
grank test indicates with high significance (P=0.96) that the two 
samples have the same properties, so that effects due to local 
compression of the ICM and source size, as well as projection 
effects should not play any role. We can then safely consider the 
sample as a whole, as we did. 

A.2. Fp trend and cluster center 

Another issue is related to the binning process. The first bin, 
t'norm < 1 could show a lowcr Fp because the number of sources 
that fall within thin bin is small, and the chance of finding 
strongly polarized sources in then reduced. To verify this, we 
repeated the same analysis described in Sec. |4]but centered on 
different position. The center is now chosen at a distance corre- 
sponding to ~10 for each cluster The Fp distribution is shown 
in Fig. lA.ll No trend is detected now. This guarantees that the 
trend detected in Sec.|4]is real and due to the cluster and not to 
missing statistics. 
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A.1. Cluster sources and background sources 

We are analyzing the Fp trend versus the cluster impact param- 
eter, so that in each bin we will have both sources that belong to 
the cluster and sources that are in its background. The difference 
in the Fp experienced by a source lying in a plane perpendicular 
to our line of sight and passing through the cluster center and a 
source that is in the background of the cluster is a small factor 
(the dispersion of the RM changes of a factor V2), so that we 
are confident that this should not cause a major effect. 
Radio sources that belong to the clusters in our sample could 
have larger angular sizes with respect to distant background 
ones. In particular, cluster radio sources sometimes show a 
"narrow-angle tail" or "wide-angle tail" morphology. Tails are 
known to be usually less polarized intrinsically, i.e. when ob- 
served at high frequencies. In addition, the ICM in the imme- 
diate surrounding of cluster radio sources could be locally com- 
pressed, so that cluster radio sources could suffer higher depolar- 
ization respect to background ones. Given that redshifts are not 
known for all of the sources in our sample, a possible way to dis- 
tinguish background and cluster radio sources is to divide them 
according to their angular size. We have separated the source 
sample into two samples having size >5 beams and < Sbeams, 
and compared the Fp trends of the two populations. The value 
of 5 beams was chosen on the basis of the radiosources in the 
cluster Abell 2255. In this cases in f act the radio sou rces that 
are cluster members are known (Mil ler & Owenll2003l) . its red- 
shift is representative of the sample, and "narrow-angle tail" and 
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Fig. A.l. Fractional polarization of sources versus projected dis- 
tance. Sources are centered on a position that is lOr^ far from 
the cluster center Arrows indicate upper limits. Bottom panel: 
crosses refer to the median of the KM estimator is each bin, bars 
indicate the 14th and 84th percentile. No trend is detected in this 
case 
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Table 1. Cluster sample selected from the HIFLUGCS catalogue 



Cluster 


RA 


DEC 


z 


Angular to linear conversion scale 


name 


(J2000) 


(J2000) 




"/kpc 


2A0335 


03h38m35.3 


+09d57m55s 


0.0349 


0.6950 


A0085 


00h41m37.8 


-09d20m33s 


0.0556 


1.0800 


A0119 


00h56m21.4 


-01dl5m47s 


0.0440 


0.8660 


A0133 


01h02m39.0 


-21d57ml5s 


0.0569 


1.1030 


A0399 


02h57m56.4 


+13d00m59s 


0.0715 


1.3630 


A0401 


02h58ni56.9 


+13d34m56s 


0.0748 


1.4240 


A0478 


04hl3m20.7 


+10d28m35s 


0.0900 


1.6790 


A0496 


03h51m58.6 


-22dl0m05s 


0.0328 


0.6550 


A0754 


09h08m50.1 


-09d38mI2s 


0.0528 


1.0290 


A 1644 


12h57ml4.8 


-17d21ml3s 


0.0474 


0.9300 


A 1650 


12h58m46.2 


-01d45mlls 


0.0845 


1.5870 


A1651 


12h59m22.9 


-04dllml0s 


0.0860 


1.6120 


A1736 


13h26ni52.1 


-27d06m33s 


0.0461 


0.9050 


A1795 


13h49m00.5 


+26d35m07s 


0.0616 


1.1880 


A2029 


15hl0ni56.0 


+05d44m41s 


0.0767 


1.4530 


A2065 


15h22m42.6 


+27d43m21s 


0.0721 


1.3730 


A2142 


15h58ml6.1 


+27dl3m29s 


0.0899 


1.6770 


A2147 


16h02ml7.2 


+ 15d53m43s 


0.0351 


0.6980 


A2163 


16hl5m34.1 


-06d07m26s 


0.2010 


3.3130 


A2199 


16h28m38.5 


+39d33m06s 


0.0302 


0.6050 


A2204 


16h32m45.7 


+05d34m43s 


0.1523 


2.6480 


A2244 


17h02m44.0 


+34d02m48s 


0.0970 


1.7950 


A2255 


17hl2m31.0 


+64d05m33s 


0.0800 


1.5100 


A2256 


17h03m43.5 


+78d43m03s 


0.0601 


1.1610 


A2597 


23h25ml8.0 


-12d06m30s 


0.0852 


1.5980 


A3558 


13h27m54.8 


-31d29m32s 


0.0480 


0.9410 


A3562 


13h33m31.8 


-31d40m23s 


0.0499 


0.9760 


A3571 


13h47m28.9 


-32d51m57s 


0.0397 


0.7860 


A4059 


23h56m40.7 


-34d40ml8s 


0.0460 


0.9040 


COMA 


12h59m48.7 


+27d58m50s 


0.0232 


0.4680 


HYDRA-A 


09hl8m30.3 


-12dl5m40s 


0.0538 


1.0470 


MKW3S 


15h21m50.0 


+07d42m32s 


0.0450 


0.8850 


ZwC11215 


12hl7m41.4 


+03d39m32s 


0.0750 


1.4240 



Col. 1 : Cluster name; Col. 2, Col. 3: Cluster center (RA, DEC); Col. 4 : redshift; Col. 5: angular scale; 
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Cluster 


L,- [0.1 -2.4keV] 




name 


lO'*" erg/s 


kpc 


2A0335 


2.53 


24!« 


A0085 


5.28 


60i 


A0119 


1.79 


365!]^ 


AO 133 


1.59 




A0399 


3.87 


332+^« 


A0401 


6.90 


182!^ 


A0478 


9.86 


73+' 

-A 


A0496 


2.03 


21+" 


A0754 


2.15 


175+'i? 


A 1644 


2.07 


218-67 
209!| 


A1650 


4.05 


A1651 


4.44 


134-f 


A1736 


1.72 


970 + 129 


A 1795 


5.49 


57! 


A2029 


9.53 




A2065 
A2142 


3.05 
11.89 


509+^^^ 


A2147 


1.54 


1-72+74 
405+1 

ioo!j 


A2163 
A2199 


20.87 
2.19 


A2204 


15.84 


5l!i 


A2244 


4.75 


93!? 


A2255 


3.04 


440+1 


A2256 


5.05 


432+28 


A2597 


3.82 


42+1 


A3558 


3.54 


163+3 


A3562 


1.67 


71+3 


A3571 


4.32 


131!^ 


A4059 


1.54 


65+3 


COMA 


4.13 


247!]^ 


HYDRA-A 


3.19 


36! 


MKW3S 


1.53 


48+j 


ZwC11215 


2.88 


319+-" 



no yS T Radio State 

lO-^cm-' keV 



c 90+0.06 A C7C0.004 ^ni+''o'' rr" 

-QQ6 0.003 -^'^ -0.07 

7 43+0.09 A C39O.OO4 f. 1 n+0.20 ppa 

-^-0.09 "•-^•'^0.004 " -0.20 

0.14!""! 0.675"°p 5.80!«| 

2.65!«;i 0.530«;S 3.80!^| CC" 

0.21+«o^ 0.713!J»'f 7.40 +«™ H" MP 

0.56!«;j;^ 0.613°flo 8.30!?i H"^ M^ 

3 26+"''' 613°°'"' 7 10+''*'' CC 

-Q.13 Q.QQ4 -Q.40 

3 q4+0.M r. 404O.O03 4 1 3+0.88 ppa 

Q 42+0.02 Q gogo.oi? g 00+'' h''3 M° 

-QQ2 "•"'°0.024 -0.05 

Q27+o.Og 579"''' 4 70+''''' M"" 

-QQ6 0024 -!!-7Q 

40+''"'* 704" '3' 5 60+''*'' D 1 

-096 flQSl -Qi9 

lOr™ 643""" 6 30+"^'' CC 

-QQ6 "•"^-'0.013 -Q.5!) 

12+'' "* 542'' '*' 3 50+'' *'' M"* 

-BQ? flQ93 -!!4Q 

2.71!; 0.596; 6.00!; CC 

3.62!»;« 0.582»;S 8.70!»;| CC 

0.19!«;JJJ 1.162°;?37 5.40!°;3° M= 

1.48!°;"^ 0.591";™ 8.80!";«" CC 

16+™ 444™ 4Ql+"-28 M' 

"•^"-0.03 " 0.046 ^-'^-0.28 ™ 

Q 44+0.02 796" °^" 13 29+" '-' H'"^ M'^' 

-90Z "•'^"0.Q28 -0.,64 " ™ 

0.81!«;i 0.655«;»^f 4.10!!J;«1 CC 

4.40!g|» 0.597°Z 7.21!;-5^ CC" 

1 QQ+0.05 607"""' 7 10+'™ CC 

-Q.Q5 Q.QtS '-^"-2.20 

0.17!«;| 0.797«;i 6.87!";^^ H^* M" 

0.25!»gl 0.914°-^^ 7.50!g:§ M' 

3 34+'':"7 633"""** 3 60+"-^" CC 

-QO6 0.008 -^•*^"-0.20 

0.44!«;«^ 0.580";""^ 5.50!!J;3;j M^ 

0.55!^;^^ 0.472;| 5.16!»!^ H^' M™ 

J Q5+0.09 613™ 6 90+'' 3'' CC 

-BQS fiQlO -03!! 

1 n+oos A C090.O01 4. io+''-3o rr" 

-Q.QS QOOi -Q.3Q 

0.29+^-^^ 0.654^™ 8.38+^1 H"** M" 

Q 4.0+0:35 r> C7^0:003 T 00+0:20 ppa 

-033 Q.QQ3 -Q.2Q 

1 91+0.24 n^si 0.008 T iTfi+O.JO r^r^a 

„,4 U.3Si„„„, ^•-'^-0.20 

0.25!-! 0.819- 5.58!- M'- 



Col. 1: Cluster name; Col. 2: X-iay luminosity in the 0. 1- 2.4 keV band; Col. 3: Cluster core radius; Col. 4: central density; 

Col. 5: jS parameter; Col. 6: cluster mean temperature; Col. 7: Radio emission. H^giant radio halo detected; Col. 8: Dynamical State. 

CC=cool-core cluster, M=merging system, D- debated, see explanation below. 

If not specified, data are from Chen et al. 2006. Quantities have been coirected for the cosmological model assumed in this paper. 



^IChen et al. (2007), 
''Iwhitakeretal. (2003 



Reiorich et al. ( 2004 
|Bpzkurt et al. ( 2009" 

Bpurdin & Mazzo ttal (l2008h. 

' Krempec-Krvgier et alj 12003) & lDicke\i fT997l) . 

= Rossetti et al. ( 200^, 

Rines&Diaferio(2006), 
' Miller & Owen (2003), 
' ISunetal. (2002), 



' iFino guenov et al. ^ 20041) . 
"Weumannetal. (2003), 
"parkevitch et al. (2003) 
pHourdin & Mazzotta (20M), 

1 no cool core according to Chen et al. (2006), but the cluster is found to be in a relaxed state with a 
moderate coo l core, by the analysis of ,Takahashi & Yamashita (2003.) 



Murgia et al. ( 2010b 
Bacchiet al. (2003 
Harris et al. ( 1980 

Kassim et al. 1 2001JX 

''^ Herbig & Birkinshaw ( 1994Y lFeretti et al.l ( 120011) 
'■^iBridle & Fomalont U976) 
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Hanisch & EricksonI ( fT98Ch 



